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In this letter, we propose a mechanism to generate large-scale magnetic fields with 
correlation lengths of 100 kpc. Domain walls with QCD scale internal structure 
form and coalesce obtaining Hubble scale correlations and align nucleon spins. 
Due to strong CP violation, nucleons in these walls have anomalous electric and 
magnetic dipole moments and thus the walls are ferromagnetic. This induces 
electromagnetic fields with Hubble size correlations. The same CP violation also 
induces a maximal hclicity (Chern-Simons) correlated through the Hubble volume 
which supports an "inverse cascade" allowing the initial correlations to grow to 
100 kpc today. We estimate the generated electromagnetic fields in terms of the 
QCD parameters and discuss the effects of the resulting fields. 



1 Introduction 

The source of cosmic magnetic fields with large scale correlations has remained 
somewhat of a mystery El. There are two possible origins for these fields: pri- 
mordial sources and galactic sources. Primordial fields are produced in the 
earlier universe, then evolve, and are thought to provide seeds which gravita- 

he fields as well 
, however, most 



tional dynamos later amplify. Galactic sources would prod 
as amplify them. Many mechanisms have been proposedB 
fail to convincingly generate fields with large enough correlation lengths to 
match the observed microgauss fields with ~ 100 kpc correlations. We present 
here a mechanism which, although probably requiring a dynamo to produce 
microgauss fields, generates fields with hundred kiloparsec correlations. We 
present this mechanism as an application of our recent understanding of QCD 
domain walls, which will be described in detail elsewhere a. 

1. Sometime near the QCD phase transition, Tqcd ~ 1 GeV, QCD domain 
walls form. 

2. These domain walls rapidly coalesce until there remains, on average, one 
domain wall per Hubble volume with Hubble scale correlations. 

3. Baryons interact with the domain walls and align their spins along the 
domain walls. 

4. The magnetic and electric dipole moments of the baryons induce helical 
magnetic fields correlated with the domain wall. 

5. The domain walls decay, leaving a magnetic field. 
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6. As the universe expands, an "inverse cascade" mechanism transfers en- 
ergy from small to large scale modes, effectively increasing the resulting 
correlation length of the observed large scale fields. 

We shall start by discussing the "inverse cascade" mechanism because it seems 
to be the most efficient mechanism for increasing the correlation length of mag- 
netic turbulence. After presenting some estimates to show that this mechanism 
can indeed generate fields of the observed scales, we shall discuss the domain 
wall mechanism for generating the initial fields and some relevant astrophysical 
phenomena associated with this mechanism. 

2 Evolution of Magnetic Fields 

As suggested by Cornwall I, discussed by Son@ and confirmed by Field and 
CarrollEl, energy in magnetic fields can undergo an apparent "inverse cascade" 
and be transfered from high frequency modes to low frequency modes, thus in- 
creasing the overall correlation length of the field faster than the naive scaling 
by the universe's scale parameter R(T). There are two important conditions: 
turbulence must be supported as indicated by a large Reynolds number Re, and 
magnetic helicity (Abelian Chern-Simons number) H = J A • Bd 3 x is approx- 
imately conserved. The importance of helicity was originally demonstrated by 
Pouquet and collaborators 13. The mechanism is thus: the small scale modes 
dissipate, but the conservation of helicity requires that the helicity be trans- 
fered to larger scale modes. Some energy is transfered along with the helicity 
and hence energy is transported from the small to large scale modes. This is 
the "inverse cascade" fflffi 

In the early universe, Re is very large and supports turbulence. This drops 
to Re ~ 1 at the e + e _ annihilation cpochp To « 100 eV. After this point (and 
throughout the matter dominated phase) we assume that the fields are "frozen 
in" and that the correlation length expands as R while the field strength decays 
as i?~ 2 . Note that the "inverse cascade" is only supported during the radiation 
dominated phase of the universe. 

Under the assumption that the. field is maximally helical, these conditions 
imply the following relationships between the initial field -B rms (Tj) with ini- 
tial correlation Z(Tj) and present fields today (T now ps 2 x 10~ 4 eV) -B rms (T now ) 
with correlation l(T now ): 

Brms(T now ) = (jT~j VTo") B rms(Ti) (1) 

/ rri \ / rji \ 5/3 



As pointed out by Sonu, the only way to generate turbulence is either by 
a phase transition Ti or by gravitational instabilities. We consider the former 
source. As we shall show, our mechanism generates Hubble size correlations 
li at a phase transition Tj. In the radiation dominated epoch, the Hubble size 
scales as T~ 2 . Combining this with (0), we see that l now cx T t 1/3 ; thus, the 
earlier the phase transition, the smaller the possible correlations. 

The last phase transition is the QCD transition, Ti = Tqcd ~ 0.2 GeV 
with Hubble size / (Tqcd) ~ 30 km. We calculate (p3) the initial magnetic field 
strength to be S rms (Ti) » eA 2 QCD /(£A QC D) ~ (IO^GV^Aqcd) where £ is a 
correlation length that depends on the dynamics of the system as discussed 
below and Aqcd ~ 0.2 GeV. With these estimates, we see that we can achieve 

B ims ~ — , I ~ 100 kpc (3) 

today. One might consider the electroweak transition which might produce 
100 pc correlations today, but this presupposes a mechanism for generating 
fields with Hubble scale correlations. Such a mechanism does not appear to be 
possible in the Standard Model. Instead, the fields produced are correlated at 
the scale T~ l which can produce only ~ 1 km correlations today. 

These are crude estimates, and galactic dynamos likely amplify these fields. 
The important point is that we can generate easily the 100 kpc correlations 
observed today provided that the fields were initially of Hubble size correlation. 
Unless another mechanism for amplifying the correlations of magnetic fields is 
discovered, we suggest that, in order to obtain microgauss fields with 100 
kpc correlation lengths, helical fields must be generated with Hubble scale 
correlations near or slightly after the QCD phase transition Tqcd- The same 
conclusion regarding the relevance of the-, QCD scale for this problem was 
also reached by Son, Field and Carroll ofl. The rest of this work presents 
a mechanism that can provide the desired Hubble size fields, justifying the 
estimate (||). We shall explain the mechanism and give simple estimates here, 
but present details of the calculations elsewhereH. 

3 Domain Walls 

The key players in our mechanism are domain walls formed at the QCD phase 
transition that possess an internal structure with QCD scak. We shall present 
a full exposition about these types of walls in another paperQ but, to be specific, 
here we shall discuss the so-called axion-77' (fly) domain wallB. 

We start with a similar effective Lagrangian to that used by Huang and 
Sikivie except that we included the effects of the rj singlet field which they 
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neglected. The Lagrangian density is 



£cff = ^|V M | 2 + ^Tr|d M U| 2 - V(U,a) (4) 

where a is the dimensionless axion field and the matrix U = exp(ir/ + in^ '\f) 
contains the pion and rf fields (to simplify the calculations, we consider only 
the SU(2) flavor group). The rj field is not light, but as we shall see, is the 
dominant player in aligning the magnetic fields so we include it. The potential 
V is given by 

V= 1 -Tr (MUe'« + h.c.) - E cos ( IM^M) (5 ) 

2 \ N c J 

which was first introduced by Halpern and Zhitnitskyi It should be realized 
that iln[det(U)] = zln(det(U)) + 2nn is a multivalued function and we must 
choose the minimum valued branch. Details about this potential are discussed 
in the original paper □ but several points will be made here. All dimension- 
ful parameters are expressed in terms of the QCD chiral and gluon vacuum 
condensates and are well known numerically: M = — diag(m* | {(fq 1 ) |) and 
E=(ba s /(32n)G 2 ). " 

The result is that two different types of axion domain walls. formEl One 
is almost identical to the one discussed by Huang and Sikivie u with small 
corrections due to the rj field. We shall call this the axion/pion (a„) domain 
wall. The second type, which we shall call the axion/eta' (a,/) domain wall 
is a new solution characterized by a transition in both the axion and rj fields 
(see our other paper □ for a complete description of this wall). The boundary 
conditions (vacuum states) for this wall are a(— 00) = rf(— 00) = and a(oo) = 
7/(00) = ±7r with 7T° = at both boundaries. 

The main difference between the structures of the two walls is that, whereas 
the an domain wall has structure only on the huge scale of m" 1 , the rf tran- 
sition in the a v i has both scales, the axion scale, TO" 1 , as well as Aqq D scale. 
Therefore, the a^i domain wall has a "sandwich" structure. The reason is 
that, in the presence of the non-zero axion (9) field, the pion becomes effec- 
tively massless due to its Goldstone nature. The rj' is not so sensitive to 8 
parameter and so its mass never becomes zero. It is crucial that the walls have 
a structure of scale Aqq D : thus there is no way for the a T wall to trap nucleons 
because of the huge differece in scales but the a v i wall has exactly this struc- 
ture and can therefore efficiently align the nucleons. The QCD domain walls 
(which were also discussed incl) have the same property as a v > walls. Namely, 
they have the structure of scale Aq^ d and they can play the same role as a v i 
walls. In what follows, for more concreteness, we use cy walls. 
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The model we propose is this: Immediately after the phase transition, the 
universe is filled with domain walls on the scale of ?q"q D . As the temperature 
drops, these domain walls coalesce, resulting in an average of Jiiie large do- 
main wall per Hubble volume with Hubble scale correlations ElO. It is these 
Hubble scale a v i domain walls which align the dipole moments of the nucleons 
producing the seed fields. 

The following steps are crucial for this phenomenon: 

1. The coalescing of QCD domain wall gives the fields rf Hubble scale 
correlations. 

2. These fields interact with the nucleons producing Hubble scale correla- 
tions of nucleon spins residing in the vicinity of the domain wall. (The 
spins align perpendicular to the wall surface.) 

3. Finally, the nucleons, which carry electric and magnetic moments (due 
to strong CP violation), induce Hubble scale correlated magnetic and 
electric fields. 

4. These magnetic and electric fields eventually induce a nonzero helicity 
which has the same correlation. This helicity enables the inverse cascade. 

4 Domain Wall Properties 

We present here a method for simplifying the calculations of the bulk properties 
of domain walls. This method makes the approximation that the domain wall 
is flat and that translational and rotational symmetries are preserved in the 
plane of the wall which we take to be the x— y plane. These approximations 
are valid in the case of domain walls whose curvature is large in comparison 
to the length scale of the pertinent physics. 

Once this approximation is made, we can reformulate the problem in 1 + 1 
dimensions (z and t) and calculate the density of the desired bulk properties 
along the domain wall. To regain the full four-dimensional bulk properties, 
we must estimate the density of the particles in the x— y plane to obtain the 
appropriate density and degeneracy factors for the bulk density. Thus, the 
final results are not independent of physics in the x—y plane, but rather, these 
effects are accounted for only through the degeneracy factors. 

4-1 Alignment of Spins in the Domain Wall Background 

We proceed to demonstrate this technique by calculating the alignment of 
fermionic spins along the wall. To estimate the strengths of the fields involved, 
we consider only the rj transition because it has a similar structure in both 
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the a v > domain wall and the QCD domain walls. We take the following simple 
interaction between the rf field and the nucleons: 



lip — m^e 



ir?'(z)75 



(6) 



For our approximations, we assume that fluctuations in the nucleon fields do 
not affect the domain walls and, thus, treat the domain walls as a background 
fieloQ. The strategy is to break (j^) into two 1 + 1 dimensional components by 
setting d x = d y = and then by manipulating the system of equations that 
result to obtain an equivalent two-dimensional system. 

First, we introduce the following chiral components of the Dirac spinors 



X2 



4-_ = 



V2 v*- 



(7) 



Secondly, we assume that our system is effectively two-dimensional (an infi- 
nite domain wall lying in the x—y plane) and hence neglect the conserved x 
and y momenta. We find that the Dirac equations which follow from are 
equivalent to the coupled system 



[id + ia 3 d 3 ] ^ + =m N e <£_ 
[id - ia 3 d 3 ] = m N e + " 1 '^+ 



(8) 



in \E'_|_ and At this stage it proves convenient to rearrange these equations 
by introducing new "two-dimensional" Dirac spinors and 



(9) 



Using the definitions (^) , we put system (JsJ) into the form of two two-dimensional 
(2D) Dirac equations 



ilndy, - m N e 171 75 



(10) 



where we have introduced the 2D Dirac matrices 

7o = o-i, 7i = ~ i<J 2, 75 = 7o7i = <?3, %% = 9^ + z»u%- 



(11) 

a A full account would take into account the effects of this back-reaction. We expect that 
they would affect the potential (pi) by altering the form of the last term _Ecos() and possibly 
adding higher order corrections, but that they would not alter the nature of the domain walls. 
Quantitatively this would alter the numerical results, but would not change the qualitative 
picture presented here. 
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Equations ( |Io| ) are reproduced from the following effective 2D Lagrangian 



C 



2D 



- m N e lv 75 



$(!) + f(2) 



i7M^M - m N e-^'fA *< 2 ), (12) 



where rj' = rj'(z) is the background classical field with boundary conditions 
rf{— oo) = 0, rj'(po) = it. This 2D Lagrangian describes two species of 2D 
Dirac fermions of oposite chiral charge intracting with the external rf field. 

Now we are ready to demonstrate that the domain walls align the spins of 
the fermions. The relevant operator (which becomes the spin operator for the 
nonrelativistic system) is 



a 
a 



(13) 



and our goal is to demonstrate that the mean value (E z ) of this operator is 
generally non-zero in the domain wall. Thus, the nucleon spins are aligned in 
the z direction and have a correlation length similar to the domain wall. 

The easiest way to demonstrate this phenomenon in our model (which is 
effectively 2D) is to use the Goldstone-Wilczek adiabatic approximation HLj 
together with a bosonization trick 



^- M sin(2Vi0i) . . 



mjv is a scale parameter). The Lagrangian ( |l2| ) after bosonization is 
1 



^2D 



-fim N [cOs(2i/7T01 - T)') + COs(2y/TT(j)2 + ??')] (15) 

The adiabatic approximation is to neglect the kinetic terms in the analysis 
of the dynamics of <f>x and </>2 fields in Equation (|l5j). In this case, the mean 
values (4>\(z)) and ($2(2)) will follow the background field r)'(z), and can be 
found by minimizing the potential U in Equation (|l5|). 



i(*)> = 



2V^ : 



>(*)> = - 



j/(f) 
20? 



(16) 



To calculate the induced spin ^^S^ in our theory we should present the spin 
operator in terms of 2D fields <j>\ and (f>2, and replace these fields by their mean 



values (16) in the domain wall background: 



-^(d z Mz)-d z Mz)), (17) 
V7T 
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where, in the last step, we used the bosonic representation for the 2D vffW 
fields^. The last step is to replace these fields by their mean- values ( ft6| ) 

($t E ^ =JVx iM!l ; (is) 

7T OZ 

where N is the appropriate normalization and degeneracy factor for the \& field 
which has canonical dimension 3/2 in four dimensions while the 2D fields 
have canonical dimension 1/2. 



4-2 Fermion Degeneracy in the Domain wall Background 

We have assumed that locally the domain walls have only a spatial z de- 
pendence. This implies that there is still a 2-dimensional translational and 
rotational symmetry in the x-y plane. These translational degrees of freedom 
imply that momentum in the plane is conserved and hence we can treat the 
neglected degrees of freedom for the fermions as free degrees. The degeneracy 
in a region of area S will simply be a sum over these degrees with a discrete 
factor g = 4 = 2 x 2 for spin and isospin degeneracy 

f dxdydp x dp y gp 2 F A qcd „ , 1(V . 

N = 9 J (27r)=» =-^r s ^— s (19) 

where we estimated the Fermi energy pp ~ h-QCD — 150 MeV. As expected, 
the degeneracy is proprtional to the area of the domain wall S. Now it is clear 
that an appropriate normalization for the two dimensional tlK*' fields can be 
achieved by adding a factor l/VS in the definition (fj). In this case these 2D 
fields have correct canonical dimension 1/2. Now we are ready to estimate the 
original four-dimensional expectation value (|l8|): 

^ '«»»M.S»«ffi.^!!£, (20) 

b ir oz ir^ oz 7r z 

which has correct dimension 3. Using the same technique one can estimate 
other matrix elements in the domain wall background which have non-zero 
magnitude and thus demonstrate that they have a large correlation L on the 
of the size of the domain wall in the x-y direction. In particular, the result for 
the mean value (^j50- xy ^) is: 

A 2 

(^ 75 *)4D ~ fJ-—^——, (21) 

7T 



6 The 2D problem under discussion is quite familiar to physics community: namely, the 
calculation of the induced fermion charge in a solitonic background. 
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where the factor ~ ^qcd/ 77 h as the same origin as in Equation ( pp[ ) and is 
related to the degeneracy of the system (|l9|), while the factor fj, ~ itln is a 
dimensional parameter originating from the bosonic representation (|l5|) of the 
effective two-dimensional theory. 

5 Magnetic Field Generation Mechanism 

Here we estimate the strengths of the induced fields in terms of the QCD 
parameters. We consider two types of interactions. First, the nucleon spins 
align with the domain wall. We assume that the fluctuations in the nucleon 
field iff are rapid and that these effects cancel, leaving the classical domain 
wall background unaltered. Thus, we are able to estimate many mean values 
correlated on a large scale on the domain walls such as (^75 o^ 1 ]/) (^IJ) and 
(^7 2 75^) ( p0| ) through these interactions as described above. These mean 
values are only nonzero within a distance Aqq D of the domain wall and are 
correlated on the same Hubble scale as the domain wall. 

From now on we treat the expectation value (^l|) as a background classical 
field correlated on the Hubble scale. Once these sources are known, one could 
calculate the generated electromagnetic field by solving Maxwell's equations 
with the interaction 

£int = -(e^cr^Ts* + (ivVio-^tyF^ + *(iD) 2 * (22) 

where (ju*) is effective electric (magnetic) dipole moments of the field iff. 
Due to the CP violation (nonzero 6) along the axion domain wall, the anoma- 
lous nucleon dipole moment in (|2|) ~ ~ is also nonzeroGJ. This is 
an important point: if no anomalous moments were induced, then only charged 
particles could generate the magnetic-field: the walls would be diamagnetic not 
ferromagnetic as argued by VoloshinEJj and Landau levels would exactly cancel 
the field generated by the dipoles. 

Solving the complete set of Maxwell's equations, however, is extremely 
difficult. Instead, we use simple dimensional arguments. For a small planar 
region of area £ 2 filled with aligned dipoles with constant density, we know 
that the net magnetic field is proportional to £ _1 since the dipole fields tend to 
cancel, thus for a flat section of our domain wall, the field would be suppressed 
by a factor of (£Aqcd) -1 - For a perfectly flai„ infinite domain wall (£ — > oo), 
there would be no net field as pointed out E_3. However, our domain walls 
are far from flat. Indeed, they have many wiggles and high frequency modes, 
thus, the size of the flat regions where the fields are suppressed is governed by 
a correlation £ which describes the curvature of the wall. Thus, the average 
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electric and magnetic fields produced by the domain wall are of the order 

- [ d *(* CT ^75*) + fi*(*«V*}] (23) 

where £* is an effective correlation length related to the size of the dominant 
high frequency modes. 

To estimate what effective scale £* has, however, requires an understanding 
of the dynamics of the domain walls. Initially, the domain walls are correlated 
with a scale of Aq^ d . As the temperature cools, the walls smooth out and 
the lower bound £(t) for the scale of the walls correlations increases from 
£(0) ~ Aqq D . This increase is a dynamical feature, however, and is thus slow. 
In addition, the walls coalesce and become correlated on the Hubble scale 
generating large scale correlations. Thus the wall has correlations from up 
to the upper limit set by the Hubble scale. Thus, the effective £* -C Hubble size 
at the time that the fields apa aligned and so the suppression is not nearly as 
great as implied by VoloshinEJ. Note that, even though the effects are confined 
to the region close to the wall, the domain walls are moving and twisted so 
that the effects occur throughout the entire Hubble volume. 

The picture is thus that fields of strength 

(E.) * (B z ) ~ -J— JLII^S™ „ £^cb 

? AqCD m N 7T £*7T 

are generated with short correlations £*, but then domains are correlated on 
a large scale by the Hubble scale modes of the coalescing domain walls. Thus, 
strong turbulence is generated with correlations that run from Aqcd up to the 
Hubble scale. 

Finally, we note that this turbulence should be highly helical. This helic- 
ity arises from the fact that both electric and magnetic fields are correlated 
together along the entire domain wall, (E) ~ (A) /r where (A) is the vector 
potential and r is a relevant timescale for the electrical field to be screened 
(we expect r ~ Aq^ d as we discuss below). The magnetic helicity density is 
thus: 

h ~ A • B ~ r(E z ){B z ) ~ r4%^- (25) 

Note carefully what happens here: The total helicity was zero in the quark- 
gluon-plasma phase and remains zero in the whole universe, but the helicity 
is separated so that in one Hubble volume, the helicity has the same sign. 
The reason for this is that, as the domain walls coalesce, initial perturbations 
cause either a soliton or an antisoliton to dominate and fill the Hubble vol- 
ume. In the neighboring space, there will be other solitons and antisolitons so 



10 



that there is an equal number of both, but they are separated and this spatial 
separation prevents them from annihilating. This is similar to how a particle 
and anti-particle may be created and then separated so they do not annihilate. 
In any case, the helicity is a pseudoscalar and thus maintains a constant sign 
everywhere along the domain wall: thus, the entire Hubble volume is filled 
with helicity of the same sign. This is the origin of the Hubble scale correla- 
tions in the helicity and in B 2 . The correlation parameter £ which affects the 
magnitude of the fields plays no role in disturbing this correlation. 

As we mentioned, eventually, the electric field will be screened. The 
timescale for this is set by the plasma frequency for the electrons (protons 
will screen much more slowly) u p which turns out to be numerically close to 
Aqcd near the QCD phase transition. The nucleons, however, also align on 
a similar timescale Aq£ D , and the helicity is generated on this scale too, so 
the electric screening will not qualitatively affect the mechanism. Finally, we 
note that the turbulence requires a^-seed which remains in a local region for 
a timescale set by the conductivity t3 a ~ cT/e 2 ~ Aqcd where for T = 100 
MeV, c w 0.07 and is smaller for higher T. Thus, even if the domain walls 
move at close to the speed of light (due to vibrations), there is still enough 
time to generate turbulence. 

For this mechanism to work and not violate current observations, it seems 
that the domain walls must eventually decay. Several mechanisms have been 
discussed for the decay of axion domain wallsli30 and the timescales for these 
decays are much larger than Aq CD , ie. long enough to generate these fields 
but short enough to avoid cosmological problems. In addition, we have found 
some additional-structures which may help solve this problem. We shall present 
these elsewhereO. In any case, we assume that some mechanism exists to resolve 
the domain wall problem in an appropriate timescale. Thus, all the relevant 
timescales are of the order Aq CD except for the lifetime of the walls which is 
substantially longer and thus, although the discussed interactions will affect 
the qualitative results, they will not affect the mechanism or substantially 
change the order of the effects. 



6 Conclusion. 

We have shown that this mechanism can generate the magnetic fields (H) with 
large correlations. It seems that galactic dynamos should still play an impor- 
tant amplification role. It seems that the crucial conditions for the dynamo to 
take place are fields B > 10~ 20 G with large (100 kpc) correlations. From ([|) 
we see that we have a huge interval 10 -10 <C £*Aqcd < 1 of £* to seed these 
dynamos. Also, if £* is small, then this mechanism may generate measurable 
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extra-galactic fields. 

We mention two new points that distinguish this mechanism from previous 
proposals First, the key nucleon is the neutron which generates the fields 
due to an anomalous dipole moment induced by the CP violating domain walls. 
The,nucleons thus make the wall ferromagnetic, not diamagnetic as discussed 
in ta. Second, the interaction between the domain walls and nucleons are 
substantial because of the similar scale (Aq^ d ) of the 77' transition in the 
domain wall. There is no way that axion domain walls with scales ~ m" 1 can 
efficiently align nucleons at a temperature Tqcd- 

The presence of the magnetic fields generated by our mechanism may have 
several observable-effects. First, large magnetic fields may alter nucleosynthesis 
production ratiostj. Secondly, lar.ge scale magnetic fields may distort the CMB 
spectrum in a measurable mannerEx These place upper bounds on the strength 
of the fields. Even the maximal fields ( p^ ) with £*Aqcd = 1 generated by 
domain walls lie within these bounds. Also, if £* turns out to be quite small, 
then, unless the distribution of galaxies is correlated with the domain walls, 
this mechanism might generate measurable extra-galactic fields. 

Two other effects may be closely related to magnetic fields generated from 
domain walls. One is the observation of ultra-high energy cosmic rays past 
the GZK cutoffEil. Magnetic fields on the scale of those discussed here may 
hold a key to explaining this mystery. The other is an apparent anisotropy 
of radiation propagation over large distances resulting in a constant offset in 
Faraday measurementsE3. One possible explanation involves the introduction 
of a Chern-Simons term by hando. This type of term might arise naturally 
from CP violating domain walls. 

Domain walls at the QCD phase transition provide a nice method of gen- 
erating magnetic fields on 100 kpc correlations today (||). In addition, the 
fields and domain walls key to this mechanism may play a role in a number of 
unexplained astrophysical phenomena. We conclude on this optimistic note. 
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